Unconventional quasiparticle lifetime in undoped graphene 
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We address the question of how small can the quasiparticle decay rate be at low energies in 
undoped graphene, where kinematical constraints are known to prevent the decay into particle-hole 
excitations. For this purpose, we study the renormalization of the phonon dispersion by many-body 
effects, which turns out to be very strong in the case of the out-of-plane phonons at the K point 
of the spectrum. We show that these evolve into a branch of very soft modes that provide the 
relevant channel for quasiparticle decay, at energies below the scale of the optical phonon modes. 
In this regime, we find that the decay rate is proportional to the cube of the quasiparticle energy. 
This implies that a crossover should be observed in transport properties from the linear dependence 
characteristic of the high-energy regime to the much slower decay rate due to the soft phonon modes. 
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The recent fabrication of single atomic layers of car- 
bon has attracted a lot of attention, as this material 
(so-called graphene) provides the experimental realiza- 
tion of a system where the low-energy electronic excita- 
tions behave as massless Dirac fermionsp], H|- The main 
results reported in Refs. Q] and can be understood 
as a reflection of the linear dependence on momentum 
of the quasiparticle energy e(k). The honeycomb lat- 
tice structure of graphene is actually known to lead to a 
bandstructure with conical shape e(k) = ±t>p-|k| around 
the corners of the hexagonal Brillouin zone, with con- 
duction and valence bands meeting at the Dirac points. 
The relativistic-like invariance arising from the massless 
Dirac quasiparticles has been shown to be at the ori- 
gin of a number of remarkable electronic properties, like 
the finite lower bound of the conductivity at the charge 
neutrality point [!, 0, IE lf|> t ne anomalous integer Hall 

effect @, HH, 

and the absence of backscattering in the 
presence of long-range scatterers@- 

More recently, the many-body properties of the 
graphene layer have been also investigated. One of the 
relevant issues addressed is whether the quasiparticle 
properties have to correspond to the expected behavior 
for a Fermi liquid in two dimensions [ToL ITl|. In this re- 
gard, it has been pointed out that the e-e interactions 
lead to quite different quasiparticle features in graphene 
depending on whether the material is doped or not [ill]. 
This can be understood from the particular kinemati- 
cal constraints of the conical dispersion, that prevent 
the Dirac quasiparticles from decaying into interband 
particle- hole excitations [12l]. 

In doped graphene, intraband processes are responsi- 
ble for the quasiparticle decay, leading to a quadratic de- 
pendence on energy of the decay rate 11]. On the other 
hand, when the Fermi level is at the charge neutrality 
point, the electron self-energy has a linear dependence 
on frequency. However, the maximum energy released in 
the scattering of a quasiparticle with momentum transfer 
q is at the boundary of the continuum of particle-hole 
excitations, which have energy > Uj?|q|. In situations 
where the Coulomb interaction remains singular in the 
limit q — > 0, as it happens in the layers of bulk graphite, 



a finite spread in the momentum of the quasiparticles is 
enough to give rise to a finite decay rate [13]]. This mecha- 
nism must rely however on some effect extrinsic to the 2D 
system (disorder, for instance), and it would be absent 
anyhow as soon as the Coulomb interaction is screened 
beyond a certain distance. 

In this paper we address the question of whether 
the quasiparticle decay rate may actually vanish in a 
graphene layer with the Fermi level tuned at the charge 
neutrality point. This study is relevant as it faces the 
possibility of having an electron liquid made of extremely 
long-lived quasiparticles. Thus, we will look for many- 
body effects which may give rise to suitable gapless ex- 
citations and consequent quasiparticle decay channels in 
undoped graphene. 

It is known that, in the absence of dopin g, th e 2D sys- 
tem does not support plasmon excitations [13J. Yet the 
polarization of the electron liquid is singular at low ener- 
gies, and this may be the source of potential instabilities. 
We will see that there is actually a significant renor- 
malization of the interactions at the large momentum- 
transfer K connecting the two inequivalent Dirac points 
in graphene. At such large momentum, the singular be- 
havior of the electron polarization tends to amplify the 
effects of the elcctron-phonon interaction, which prevails 
over the Coulomb interaction. Thus, we will show that 
gapless phonon branches may appear at the K point 
when graphene is lying on a substrate. The resulting low- 
energy phonon modes provide then the relevant mecha- 
nism for the decay of quasiparticles in undoped graphene, 
though with a very low decay rate that turns out to be 
proportional to the cube of the quasiparticle energy. 

We begin by considering the hamiltonian for Dirac 
quasiparticles in graphene, at energies below the scale of 
~ 1 eV for which the dispersion can be taken as linear: 



H = v F J d 2 fc* (a)t (k) 7 < a ) -k*W(k) 



(1) 



In the above expression, a sum is implicit over the index 
a accounting for the two different valleys and correspond- 
ing Dirac spinors at opposite corners K,—K in the 
graphene Brillouin zone. 7 < - a - 1 are different sets of Pauli 
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matrices for a = 1,2, which must be chosen according 
to the appropriate chirality of the modes at K, —K as 
7W = (a x ,<T y ), 7 (2) = (-a x ,a v ) 0- As a first step in 
the development of the many-body theory, we will as- 
sume that the quasiparticles interact through a Coulomb 
potential 



2k 4 



(2) 



with a dielectric constant k dictated by the coupling to 
the substrate. 

As is well-known, the quasiparticles of the 2D system 
provide very limited screening of the Coulomb potential 
in @. This effect can be assessed by computing the 
polarization 



4Tr 



d k f dZJk 
"2V 



(2tt) 2 

G<°> (k + q, iu3 k + ioj q ) G [h] (k, iuj k \3) 



with Dirac propagators G^ a ^(k, iZJk) = — Vp^ a ' ■ 

k). At small momentum-transfer, the trace in is taken 
over excitations in the same valley a = £>, with the result 
that Tr(iZU g + vf'y'"'' ■ q)(«^fc + vpf^ ■ k) = —2UJ q uJk + 
2vpq ■ k. This leads to an expression for the intraval- 

ley polarization r±Q a ' a ^(q, iU3 q ) = — q 2 



vpq 2 



^ .1- 



Going back to real frequency uj q = icJ q , we find a diver- 
gence of the polarization at ui q = Vp\q\. This marks actu- 
ally the threshold for the creation of particle-hole pairs 
in the electron liquid. The particle-hole continuum is 
above the maximum energy vp\q\ released in the scatter- 
ing of a quasiparticle with momentum transfer q. This 
explains that the quasiparticle decay into particle-hole 
pairs is forbidden in the case of undoped graphene. 

In the case of intervalley scattering of quasiparticles, 
the polarization is also affected by a similar divergence 
at io q — i>j?|q|, where q stands now for a small devi- 
ation around the large momentum K . The computa- 
tion of the polarization ^ with a 7^ b leads to the 
trace Tr(iU q + vp^f^ ■ q)(iwfc + vpj^ ■ k) = — 2Lu q uJk — 
2vpq x k x + 2vpq y k y . This can be assimilated to the 
above computation for a = b if the y component of each 
momentum is exchanged with the frequency to, and an 
overall — sign is introduced. It can be checked by di- 
rect calculation that the result for the intervalley po- 
larization corresponds actually to operating that trans- 
formation in the above expression for IlQ a ' a \ that is, 

n<, 1,2) (q,«U g ) = (p 2 +LJ 2 q /vp)/8^v F c?+LJ 2 q . The po- 
larization thus obtained shows a preferred direction in 
momentum space, which is a reflection of having con- 
sidered the scattering between two Dirac valleys along 
the x direction. The result physically sensible can be ob- 
tained by averaging over the processes involving the three 
equivalent nearest-neighbor valleys of the K point. These 
include in particular the valleys rotated by an angle of 
±27r/3 with respect to the :r-axis. Taking into account 



the three different contributions, we get the final result 
for the intervalley polarization 



n (q,w 9 



q 2 /2-^A,f. 



(4) 



The polarization ^ does not give rise to any singu- 
larity when renormalizing the Coulomb interaction, as 
the Coulomb potential gets dressed at large momentum- 
transfer K in the form Vb(q) ~ c 2 /(2kK — Ho(q,u) q )). 
On the contrary, the singular behavior of the interval- 
ley polarization may lead to important effects in the 
phonon sector. This consideration is relevant for lat- 
tice vibrations coupling to the total electron charge, as 
it happens in the case of the out-of-plane phonons. The 
electron-phonon interaction can be analyzed in terms of 
the atomic deformation potential induced by the lattice 
vibrations [15|. When graphene is lying on a substrate, 
the mirror symmetry of the vibrations perpendicular to 
the carbon layer is broken, and the on-site deformation 
potential induces a linear coupling of the out-of-plane 
phonons to the total electron charge. If we denote by 
c q the creation and annihilation operators for the 



modes around the K point of a branch of out-of-plane 
phonons, we can describe the phonon sector by means of 
kinetic and interaction terms in the hamiltonian: 

H ph = J d 2 qui (q) c q c q 

ffc-ph = 9 J d 2 kd 2 q^(k + (l )^(k)(c a + c^)(5) 

For the characterization of the phonon branch, it will be 
enough to approximate the energy of the out-of-plane 
phonons about the K point by uiq m 70 meV. The 
electron-phonon coupling g can be obtained as the atomic 
deformation potential (of the order of a few eV) times 
1 / y/mcUJo (mc being the carbon atomic mass) [15| . 

The intervalley polarization induces a strong renormal- 
ization of the out-of-plane phonons in graphene. It has 
been already pointed out that the interaction with the 
electronic degrees of freedom may give rise to signifi- 
cant Kohn anomalies in the dispersion of in-plane op- 
tical phonons [Lsl \vi\ . The coupling to these branches is 
given in general by the modulation of the transfer inte- 
gral between nearest-neighbor atoms in the carbon lat- 
tice. In our framework, this gives rise to an electron- 
phonon vertex proportional to the matrix a x . When in- 
troduced in the computation of the trace in the polariza- 
tion, such a vertex gives rise to simple scalar products in 
(q,!Ug) space, leading to a susceptibility proportional to 

J vpq 2 — a-> 2 . However, in the case of phonons coupling 

to the total electron charge, the particle-hole polarization 
(|4]) induces a more profound anomaly in the phonon dis- 
persion. If we approximate the bare phonon propagator 
about the K point by D (q,u}) w 2Lu a /(ui 2 — uu 2 + ie), 
the renormalized propagator D(q,w) dressed with the 
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FIG. 1: Plot of the region corresponding to the particle-hole 
continuum (shaded area) and the phonon dispersion arising 
from the solution of Eq. 0. 



particlc-holc polarization becomes 

2uj 



D(<i,lj) 



2wo.g 2 n (q,w) 



(G) 



The renormalizcd phonon energies are found by setting to 



zero the denominator of the propagator 
to the equation 



which leads 



U3 



g 2 v 2 F q 2 /2 

— UJQ 

F 



V 



= o 



(7) 



It can be checked that the phonon dispersion thus com- 
puted becomes gapless at q = 0, adopting the form of a 
low-energy branch below the continuum of particle-hole 
excitations as shown in Fig. [TJ At this point, it becomes 
pertinent however to assess the effects of the Coulomb 
interaction on the renormalization of the phonon proper- 
ties. 

The Coulomb and the phonon-mediated interaction 
have to be considered on the same footing when ana- 
lyzing their role in the renormalization of the phonon 
propagator. Thus, we can define an intervalley particle- 
hole susceptibility II(q, w q ) dressed by the effect of the 
Coulomb interaction and, therefore, satisfying the equa- 
tion 



n = n n 



n i/ n 



(8) 



In terms of this susceptibility, the renormalized phonon 
propagator D(q, u>) can be obtained as 



D- 1 = Dq 1 - g 2 Tl 



(9) 



Combining the solution of 
final expression 



D(q,u) 



The strength of the phonon-mediated interaction is given 
by the dimcnsionlcss coupling g 2 /vp 7 which can be esti- 
mated as ~ 0.1. This is smaller than the strength of the 
Coulomb interaction in typical graphene samples, where 
e /kvf ~ 1. However, the latter enters in the denomina- 
tor of pop with a relative weight ujq/AvfK, which is of 
the order of ~ 0.001. Therefore, we see that the Coulomb 
interaction is not able to balance the effect of the strong 
renormalization of the phonon propagator at the large 
momentum transfer K. 

We find then that, when graphene is tuned at the 
charge neutrality point, there are phonon modes with 
very low energy around the K point, arising as a conse- 
quence of the strong renormalization from the coupling to 
particle-hole excitations. Strictly speaking, the phonon 
branches become gapless only at zero temperature in the 
undoped system. Away from the charge neutrality point 
and at finite temperature, the divergence of the polariza- 
tion in will be cut off by either the thermal energy 
or the effective chemical potential of the system. Above 
such infrared scales, the phonon dispersion will follow 
anyhow the trend represented in Fig. [TJ The phonon en- 
ergy Wph(q) obtained from the renormalized propagator 
is actually given by 



(11) 



with 10), we arrive at the 

2^o(l- ^n (q,^)) 

- + it - ({lu 2 - w§) 3^ + 2w o5 2 )n (q, u) 

(10) 



^ph(q) « v F \q\ - — T 



8v 2 F 



2w§ 



We stress that this phonon branch lies in any event away 
from the continuum of particle-hole excitations, opening 
the possibility to observe well-defined phonon modes of 
very low energy at the K point of graphene. 

The existence of the soft phonon branch leads to 
a channel for the decay of quasiparticles in undoped 
graphene. The maximum energy that can be released 
by a quasiparticle in a scattering process at any low 
momentum-transfer is enough to hit the phonon branch 
pip, so that electron quasiparticles can decay into this 
type of phonon modes down to arbitrarily low energies (at 
zero temperature). The quasiparticle decay rate r _1 can 
be computed from the electron self-energy £( Q )(k,a;fc) as 

t- 1 = -Im S (a) (k,w F |k|) 

« Im z ff 2 J ^^G< b >(k - q, v F \k\ - W ,)r>(q, Uq ) (12) 

which amounts to making the convolution of the imag- 
inary part of the electron propagator with that of 
D(q,u q ). 

In Eq. (fT2|) . the imaginary part of enforces the 
constraint uj q — Vp\k.\ — Vp\\c— q|. For that frequency, the 
phonon propagator picks up an imaginary contribution 
only from the phonon branch (fTTj) . We have actually 



d 2 g 

(27T) 2 



<J(Q(q,fi q )) 



(13) 
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where fL 



UfW — vp\]i — q| and 



Q( q , n q 



ft 2 



w 2 



5 2 V 2 q 2 /2-f2 2 



2lu 



v f 8JvW 



Q 2 



(14) 



The integral over q can be done by trading the azimuthal 
variable of integration <f> by f2 q . Thus we get 

»|k| 



1 r\K-\ /-|q| 

1 « 5-^ / ^ |q| / dO, 
2tt Jo Jo 



9Q 



5(fi q - Wph(q)) (15) 



The expression (jl5j> leads to different behaviors de- 
pending on whether the quasiparticle energy is well above 
or below the scale u)q. In the range where vf\^\ t^o, 
it is easy to see that the Jacobian |<9(/>/<9£l q | scales as 
- |k|/|q| v /4k 2 -q 2 , while \dQfdflq\- 1 does not scale 
with momentum. The quasiparticle decay rate shows 
then a linear dependence on energy r _1 ~ (g 2 /v F )vp\k\, 
in agreement with previous analyses of the decay due to 
optical phonons 18] . On the other hand, when the quasi- 
particle energy is below u)q, we find that \d<f)/dClq\ scales 

as ~ cdoVW ~ Iql/VW 11 ^ 2 - Moreover, we also have 
\dQfdQq\~ 1 ~ v F \q\ 3 /u>Q. We get then a decay rate 



1 
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167T VF 



dx x yl 



(16) 



We arrive at the result that, in the case of undoped 
graphene, the quasiparticle decay rate cannot vanish at 
low energies, even below the scale lvq of the out-of-plane 
phonons, where it must be proportional to the cube of 
the quasiparticle energy. 

The behavior (fT6|) differs significantly from the rate 
obtained for a long-range Coulomb interaction, which is 
proportional to the quasiparticle energy [Til [l3j . We re- 
mark that the two behaviors correspond actually to quite 
different conditions. When the Coulomb interaction re- 
mains long-ranged, as in the layers of bulk graphite, the 
singular character of the potential @ leads to a jump 
in the electron self-energy at uik = vf\^\- A spread in 
momentum of the quasiparticles (induced for instance by 
disorder) may be invoked to obtain a finite quasiparticle 
decay rate by taking the limit ojk — vp\k\ + + . We have 
to bear in mind however that, in graphene, the divergence 
of the Coulomb potential may be cut off by some finite 



screening length I. If we assume a potential of the form 
Vb(q) = e 2 /-\/q 2 + l~ 2 , we can estimate the quasiparti- 
cle decay rate to lowest order in the Coulomb interaction 
as 



lim e 4 



dq |q| 



q|+e 



|q| 



VMllV 6 - (ft q -UF|q|) ( q 2 + /-2) / n a _ v 2 q2 



(17) 



We see that the decay rate in this approach turns out to 
be proportional to (e 2 /vp) 2 vpl 2 \k\ 3 . It is only after send- 
ing the screening length to infinity that the divergence of 
the Coulomb potential at q = is able to change the 
scaling of the integrand in (|17[) . leading to a decay rate 
proportional to the quasiparticle energy. 

In conclusion, we have shown that the dispersion of 
phonons that couple to the total electron density under- 
goes a strong renormalization when graphene is tuned to 
the charge neutrality point. This reflects in the behavior 
of the out-of-plane phonons, which turn out to develop a 
gapless branch at the K point of the spectrum. In situa- 
tions where graphene is very lightly doped, the mentioned 
renormalization will still give rise to a branch of very soft 
phonons, with a small gap proportional to the effective 
chemical potential (as measured from the charge neutral- 
ity point). We have seen that this branch provides the 
relevant channel for the decay of quasiparticles below the 
typical scale of the out-of-plane phonons, lu ~ 70 meV. 
Thus, even at such low energies, the quasiparticle decay 
rate does not vanish in undoped graphene, though it be- 
comes very suppressed, with a dependence proportional 
to the cube of the quasiparticle energy. The present 
analysis may then be useful to interpret the results of 
transport experiments in undoped or very lightly doped 
graphene, where a crossover should be observed from the 
linear decay rate characteristic of the high-energy regime 
to the much slower cubic dependence found in the paper. 
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